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Abstract. Topological superconductors have become a subject of intense 
research due to their potential use for technical applications in device fabrication 
and quantum information. Besides fully gapped superconductors, unconventional 
superconductors with point or line nodes in their order parameter can also exhibit 
nontrivial topological characteristics. This article reviews recent progress in the 
theoretical understanding of nodal topological superconductors, with a focus 
on Weyl and noncentrosymmetric superconductors and their protected surface 
states. Using selected examples, we review the bulk topological properties of 
these systems, study different types of topological surface states, and examine 
their unusual properties. Furthermore, we survey some candidate materials for 
topological superconductivity and discuss different experimental signatures of 
topological surface states. 
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1. Introduction 

The discovery that gapped electronic systems can possess nontrivial topology has 
sparked a revolution in condensed matter physics [T|[2|5|[11|5|E1[7|[H]. Such 
topological systems are characterized by an Z or Z 2 topological number, and 
display topologically protected surface states due to the so-called bulk-boundary 
correspondence. Two species of such systems can be discerned: topological insulators 
and superconductors. While there are now many examples of topological insulators, 
topological superconductors have proved to be much rarer. This is due to the 
unconventional pairing symmetries required for a topological state. In contrast 
to the familiar s-wave spin-singlet pairing in conventional superconductors, such 
pairing states appear to require a particular constellation of pairing interactions and 
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electronic structure, and are not robust to disorder. Topological superconductors 
are nevertheless of great interest due to the existence of protected Majorana 
surface states, which arise as a consequence of a nontrivial topology of the bulk 
quasiparticle wave functions. These surface states may have important applications in 
quantum information technology, and have motivated an intense effort to “engineer” a 
topological superconducting state from more prosaic components [gilOllIIlIE]. While 
there are many such proposals, at present the most promising is to proximity-induce 
p-wave superconductivity in a semiconductor with strong spin-orbit coupling [nm 

usmniin]- 

Although superconductors with unconventional pairing symmetries are not 
uncommon, they rather typically display points or lines on the Fermi surface where 
the superconducting gap vanishes, so-called “nodes”. As such, they are not covered by 
the “ten fold way” classification of the topology of fully gapped electronic systems [U 
HI n m. Nevertheless, such systems can display interesting topological properties, 
and possess topologically-protected edge states, although the protection is weaker 
than for fully-gapped systems [iaiiniiiniiiiiiiiH3!MiH3iiaiiaHHi!2iii3niis5- 
Although global topological invariants cannot be defined for gapless systems, it is 
nevertheless possible to classify the rich topological structure of nodal superconductors 
in terms of momentum-dependent topological numbers. This has revealed that many 
nodal topological superconductors and superfluids have nontrivial topology, e.g. the 
high-Tc dx 2 -y 2 -wave superconductors [32 |331 [Ml |3S], the heavy Fermion systems 
CeCoIns |36l Ell EH El] and URu 2 Si 2 [401 lU IMl HSj, nodal noncentrosymmetric 
superconductors |2Q1 ED ES ESI ESI ESI ES] , Weyl superconductors |46l ED ESI ES] , 
oxide interfaces [^, and the A phase of superfluid ^He [12 ESI ED ESI ED ■ Similar 
to fully gapped topological superconductors, the topological characteristics of nodal 
superconducting states reveal themselves at the surface in the form of Majorana cone 
states, Majorana flat-band states, or arc surface states |22j . 

Our aim in this review is to provide a thorough introduction into the held 
of nodal topological superconductors, acting as a manual and reference for both 
theorists and experimentalists. We do not wish to deliver an exhaustive overview 
of all nodal superconductors, however, but rather to show via selected examples the 
different topological structures which can be realized and what type of unconventional 
properties theses materials exhibit. While the concepts that we discuss can be 
generally applied to all nodal topological systems, our focus shall be upon systems 
where the superconductivity is an intrinsic property, as opposed to engineered nodal 
superconducting states in heterostructures. 

We give a brief outline of the review here. We commence in Sec. E by 
introducing our two model systems: Weyl superconductors and noncentrosymmetric 
superconductors, which are prototypes for topological superconductors with point 
and line nodes, respectively. In Sec. we set out the topological classification of the 
nodal states in terms of a “periodic table” of the nodal structure, and apply this 
to our canonical examples. The topologically-protected surface states are studied 
in Sec. using a quasiclassical method, which forms a useful basis for understanding 
the experimental signatures of the topological surface states in Sec. We survey 
candidate nodal topological superconductors in Sec. before concluding with a brief 
outlook for the held. 
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2 . Phenomenological models 

While there are a plethora of nodal superconductors with nontrivial topology, it is 
convenient to illustrate our discussion with straightforward but physically relevant 
phenomenological models. Here we introduce three such systems: the p-wave 
Weyl superconductor, the d-wave Weyl superconductor, and noncentrosymmetric 
superconductors. Our choice is motivated to exemplify line and point nodal structures, 
and we shall return to these model systems throughout the review. 

S.J. Topological superconductors with point nodes 


We consider two different phenomenological models of superconductors/superfluids 
with point nodes. The first one is the Anderson-Brinkman-Morel (ABM) state |52l 
[531 ED, which is believed to occur in the A phase of superfluid ^He as well as in cold 
atomic and polar molecular gases |^. The ABM state is a three-dimensional super¬ 
fluid of spinless (or fully spin-polarized) fermions with ± ipy)-'waYe pairing and 
is described by the Bogoliubov-de Gennes (BdG) Hamiltonian TL = 5 
with 

'H(k) =/i(k)r^-f • (eiTj,-b e2Ty), (1) 


where (/)k = (ckjcj^.)"*" and c. 


( Cj. ) is the creation (annihilation) operator of spinless 
fermions with momentum k. Here, denote the three Pauli matrices which act in 
particle-hole space, 1 = ei x 62 represents the direction of the orbital momentum of the 
Cooper pair, /i(k) is the normal-state dispersion, kp is the Fermi wave vector, and Aq 
denotes the amplitude of the superconducting order parameter. In the following we 
assume a spherical Fermi surface centered at the F point and a momentum independent 
vector 1 pointing along the axis in the Brillouin zone. Thus, the eigenvalues of 
Hamiltonian (1) are ±A(k) = ±y^[/i(k)]2 -|- {Ao/kF)'^{k'^ + ky), and the spectrum 

A(k) exhibits two Weyl nodes at the north and south poles of the Fermi sphere 
K± = (0,0,±fcF)- The low-energy physics in the vicinity of the point nodes K± 
can be captured by the anisotropic Weyl Hamiltonian 


'H±{p) = {Ao/kF){pj;T^ NPyTy) ± hvpPzTz, ( 2 ) 

where p = k—K± and vp = hkp/rn denotes the Fermi velocity. Hence, Eq. Q defines 
the simplest version of a Weyl superconductor with two point nodes. The dispersion 
close to the Weyl points is linear in all three directions, leading to a density of states 
that increases quadratically with energy. As opposed to Weyl points in semimetals, 
the particle-hole symmetry dictates that the point nodes K± remain pinned at zero 
energy even as the chemical potential is varied. Each Weyl node can be characterized 
in terms of a chirality index, which measures the relative handedness of the three 
momenta p with respect to the Pauli matrices in Eq. (H). Eor Hamiltonian 'H± (p) the 
chirality of the Weyl point is ±1. We note that while arbitrary (small) perturbations 
of Hamiltonian (p|) in general lead to a shift of the point nodes, they cannot open up 
a gap in the spectrum. We will see in Sec. |3.3.l1 that the stability of the Weyl nodes is 
protected by a Chern number (or Skyrmion number), which takes on the values ±1. 

As a second example of a superconductor with point nodes, we consider the three- 
dimensional spin-singlet chiral {d^ 2 _y 2 ± idxy)-YfaYe state, which is likely realized 
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in the heavy fermion system URu 2 Si 2 [40l |4TJ |42l |43] and the pnictide material 
SrPtAs [m [37j . The BdG Hamiltonian of this chiral superconducting phase is given 


by^ = iEk‘&ic^(k)$k: 


with 


'H(k) = h(k)r^ 


kl 


ky)Tx + “^kxkyTy^ 


(3) 


and the Nambu spinor <i)k = (Ck ^^j where cl^ (Ck^) represents the electron 

creation (annihilation) operator with momentum k and spin s. Observe that for this 
pairing symmetry the orbital angular momentum projected along the kz axis is ±2, 
whereas for the {px ± ipy)-wa.ve state discussed above it is ±1. Similar to Eq. 
we find that for a spherical Fermi surface centered at the P point, Hamiltonian 
exhibits two point nodes located at K± = (0,0, T/cf), which are double Weyl points 
pinned at zero energy by the particle-hole symmetry. Expanding Eq. ([^ around Kj- 
to second order in the momentum, we obtain the following low-energy description of 
the double Weyl nodes 


'H±{P) = [{pI-pI)Tx + ^PxPyTy] + 


2m 


p| ±hvFPz 


(4) 


where p = k —Kj-. Hence, the low-energy nodal quasiparticles near Kj- exhibit linear 
and quadratic dispersions along the Pz direction and in the PxPy plane, respectively. All 
in all, this anisotropic dispersion leads to a density of states which increases linearly 
with energy. As explained in Sec. 3.3.2 the double Weyl points Q are protected 
against gap opening by the conservation of a Chern or Skyrmion number that takes 
on the values ±2. 


2.2. Topological superconductors with line nodes 


The prototype of a topological superconductor with line nodes are time-reversal- 
symmetric three-dimensional noncentrosymmetric superconductors (NCS). In these 
systems, the crystal structure breaks inversion symmetry leading to strong electric 
held gradients within the unit cell, which in turn generates strong spin-orbit coupling 
(SOC). Furthermore, the absence of inversion symmetry means that parity is no longer 
a good quantum number, and so a mixed pairing state with both singlet and triplet 
gaps is possible. 

The minimal description for such a system has the BdG Hamiltonian H = 
5Ek^k'^(k)^k: with 

/ h{k) A(k) \ 

At(k) -h'^i-k) ) 

and ipk = (Ck^) Ck^, c^kf''-k^)"'"> where Ck^. (Ck^,) denotes the electron creation 
(annihilation) operator with momentum k and spin cr. The normal-state dispersion 
of the electrons in the spin basis is given by 

h(k) = ekCTo + gk • cr, (6) 

where £k is the nonmagnetic dispersion while gk denotes the SOC potential, cr = 
{crx,cry,az) denotes the vector of the Pauli matrices, while (Tq stands for the 2x2 
unit matrix. By time-reversal-symmetry, we require that Ek and gk are symmetric 
and antisymmetric in k, respectively. The normal-state Hamiltonian Eq. § is 
diagonalized in the so-called helicity basis, taking the form /^(k) = diag(^j^, ^j^), where 
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= Ek ± |gk| is the dispersion of the positive-helicity and negative-helicity bands, 
respectively. 

The superconducting gap contains both even-parity spin-singlet and odd-parity 
spin-triplet pairing potentials 

A(k) = (V'kO'o -I- dk • cr) {i(Ty) , (7) 

where ^pk and dk represent the spin-singlet and spin-triplet components, respectively. 
In a time-reversal-symmetric system they have the same global phase, and without 
loss of generality we assume them to be real. It is well known that in the absence 
of interband pairing, the superconducting transition temperature is maximized when 
the spin-triplet pairing vector dk is aligned with the polarization vector gk of the 
SOC [55]. This implies that there is only pairing between states on the same helicity 
Fermi surface, with positive and negative helicity gaps = ^k ± |dk|. Assuming an 
s-wave like order parameter for the singlet pairing, we observe that while the positive- 
helicity Fermi surface is fully gapped, the negative-helicity gap may possess nodes. We 
parametrize the singlet and triplet components of the superconducting gap function 
as 

ipk = As = rAo, (8) 

dk = Atlk = (I - r) Aolk, (9) 

where Ik = gk/A, with A the SOC strength. Here, r parametrizes the mixed pairing 
state, tuning the system between purely triplet (r = 0) and purely singlet (r = 1) 
gaps. 

3. Topological classification of nodal superconductors 

The topological properties of nodal superconductors (and semimetals) can be classified 
in a similar manner as those of fully gapped superconductors (and insulators) [H [H [3] . 
A nodal point or nodal line in a superconductor is topologically stable if there 
does not exist any symmetry preserving mass term that opens up a full gap in 
the spectrum. The topological stability of these point or line nodes is guaranteed 
by the conservation of a topological charge (i.e., a topological invariant), which is 
defined in terms of an integral along a contour that encloses (encircles) the point node 
(line node). These topological invariants can be either Chern or winding numbers, 
referred to as “Z invariants” (or “2Z invariants”), or binary invariants, referred to 
as “Z 2 numbers”. In this section we briefly review the ten-fold classification of 
nodal topological superconductors (and semimetals) [SS] (57] [551 IHl HO] • This scheme 
classifies point and line nodes in nodal superconductors (Fermi points and Fermi 
lines in semimetals) in terms of global symmetries, i.e., nonspatial symmetries that 
act locally in position space, namely, time-reversal symmetry (TRS), particle-hole 
symmetry (PHS), and chiral or sublattice symmetry (SLS). 

In momentum space, TRS and PHS act on the BdG (or Bloch) Hamiltonian 77 (k) 
as 

T-^n{-k)T =+n{k) and C-^n{-k)C =-n{k), (10) 

respectively, where T = ICUt and C = ICUc are antiunitary operators and K. denotes 
the complex conjugation operator. Chiral symmetry, however, is implemented in terms 
of an anticommutation relation 

s-^n{k)s= -n{k), 


( 11 ) 
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where S' is a unitary operator. Chiral symmetry S can be viewed as a combination 
of time-reversal and particle-hole symmetry, i.e., S (x TC. There are three different 
possibilities for how a given BdG (or Bloch) Hamiltonian 77(14) can transform under 
TRS (or PHS): (i) 77(14) is not symmetric under TRS (or PHS), which is denoted by 
“0” in Table[^ (ii) 77(k) is invariant under TRS (or PHS) and T (or C) squares to -|-1, 
which is denoted by in Table[^ and (iii) 77(k) is symmetric under TRS (or PHS) 
and T (or C) squares to —1, which is denoted by “—1” in Table Hence, there are 
3x3 = 9 possibilities for how a Hamiltonian can transform under both TRS and PHS. 
For eight of these cases, the presence or absence of SLS is fully determined by how 77 (k) 
transforms under TRS and PHS. However, when both time-reversal and particle-hole 
symmetry are absent, there exists the extra possibility that the combined symmetry 
S 04 TC is still present. Hence, altogether there are ten cases, that fully exhaust 
all the possible transformation properties of a given Hamiltonian 77 (k) under global 
symmetries. These ten cases, each of which defines a symmetry class, are listed in the 
first column of Table We note that these ten symmetry classes are closely related 
to the ten classes of random matrices introduced by Altland and Zirnbauer HD Eg. 

Now, the ten-fold classification of nodal superconductors (and semimetals) does 
not only depend on the symmetry properties of 77 (k) but also on the codimension 

P = dez — dn (12) 

of the nodal structure, i.e., on the difference between the dimension dez of the Brillouin 
zone and the dimension of the superconducting node. Furthermore, it depends on 
how the superconducting nodes transform under the nonspatial symmetries ED] . 
Two different situations have to be distinguished: (i) Each point node or line node 
is left invariant by nonspatial symmetries and (ii) different point nodes or line nodes 
are pairwise related to each other by global symmetry operations. In case (i) the 
superconducting nodes are located at high-symmetry points of the Brillouin zone, i.e. 
at time-reversal invariant momenta, whereas in case (ii) the point or line nodes are 
positioned away from high-symmetry points of the Brillouin zone. 

3.1. Nodes at high-symmetry points 

Let us start by discussing the ten-fold classification of superconducting nodes located 
at high symmetry points of the Brillouin zone [SDl ED EHl ED ED]- l^is case 
the classification of point nodes can be derived from the ten-fold classification of 
fully gapped superconductors trough a dimensional reduction procedure [5D|- That 
is, the surface state of a fully gapped dez-dimensional topological superconductor 
realizes a topologically stable point node in dez — 1 dimensions. Indeed, the bulk 
topological invariant of the dez-dimensional fully gapped superconductor is directly 
related to the topological charge of the nodal point at the boundary EH E3] ■ From 
these observations it follows that the classification of global-symmetry-invariant nodal 
points with d„ = 0 is obtained from the original ten-fold classification of fully gapped 
superconductors by the dimensional shift dez —t dez — I- This reasoning also holds 
for nodal structures with codimension p < dez (i-e., dn > 0), provided that they are 
protected by a Z invariant or 2Z invariant. Z 2 numbers, on the other hand, guarantee 
only the stability of nodes with d^ = 0, i.e., point nodes. These findings are confirmed 
by derivations based on K theory ED ED ED and minimal Dirac Hamiltonians EOj. 
The classification of global-symmetry-invariant nodal structures (and Fermi surfaces) 
is summarized in Table [D where the first row indicates the codimension p of the 


CONTENTS 


Table 1. Ten-fold classification of nodal points/lines and Fermi surfaces in 
superconductors and semimetals, respectively |56l [ 5 TI [58l |60|. The first and 
second rows indicate the codimension p = d^z ~ of tho nodal lines at high- 
symmetry points and away from high-symmetry points of the Brillouin zone, 
respectively. The first column lists the ten symmetry classes (using the “Cartan 
nomenclature” m), which are distinguished by the presence or absence of TRS 
(T), PHS (C), and chiral symmetry (S). The second and third columns indicate 
the sign of the squared symmetry operators and , respectively. The absence 
of symmetries is denoted by “0” . The presence of chiral symmetry S is denoted 
by “1”. 
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^ Z 2 numbers only protect nodal lines of dimension zero (dn = 0) at high-symmetry 
points of the Brillouin zone. 


§ . . . 

Nodal lines of any dimension dn located away from high symmetry points of the 
Brillouin zone cannot be protected by a Z 2 number. However, the Z 2 topological 
invariant can lead to the appearance of protected zero-energy surface states at time- 
reversal invariant momenta of the surface Brillouin zone. 


superconducting nodes. For any codimension p there are three symmetry classes for 
which stable superconducting nodes (or Fermi surfaces) exist that are protected by a 
Z invariant or 2Z invariant, where the prefix “2” indicates that the topological number 
only takes on even values. Furthermore, in each spatial dimension dez there exist two 
symmetry classes that allow for topologically stable point nodes (Fermi points) which 
are protected by a binary Z 2 number. 


3.2. Nodes off high-symmetry points 

Second, we review the topological classification of superconductors with nodes that 
are located away from high-symmetry points of the Brillouin zone [5S1 ISO] . These 
point or line nodes are pairwise mapped onto each other by the global antiunitary 
symmetries, which relate k to —k. An analysis based on the minimal-Dirac- 
Hamiltonian method m shows that only Z-type invariants can guarantee the stability 
of superconducting nodes off high-symmetry points, whereas Z 2 numbers do not lead 
to stable nodes. However, as illustrated in terms of the example of Sec. 3.3.3[ Z 2 
numbers may nevertheless give rise to zero-energy surface states at time-reversal- 
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Figure 1. Energy spectrum of the three different types of topological subgap 
states that can exist at the surface of nodal noncentrosymmetric superconductors: 
(a) helical Majorana cone, (b) arc surface state, and (c) flat-band surface state. 
Figure adapted from Ref. m- 


invariant momenta of the surface Brillouin zone. The complete classification of 
superconducting nodes that are located away from high-symmetry points is presented 
in Table[^ where the second row gives the codimension p of the superconducting node. 
Observe that this classification scheme is related to the ten-fold classification of fully 
gapped superconductors and insulators by the dimensional shift dez dez + 1- 

3.3. Examples 

For the phenomenological model Hamiltonians given in Sec. we derive in 
this subsection explicit expressions for the topological invariants that protect the 
superconducting nodes against gap opening. We also use these examples to illustrate 
the bulk-boundary correspondence [531 El] j which links the topological characteristics 
of the nodal gap structure to the appearance of zero-energy states at the boundary. 
Depending on the case, these zero-energy surface states are either linearly dispersing 
Majorana cones, Majorana fiat bands, or arc surface states (cf. Fig. [^. We note that 
in real superconducting materials the gap nodes are usually positioned away from the 
high-symmetry points of the Brillouin zone. Indeed, this is the case for the three 
examples of Sec. which are therefore classified according to Sec. |3.2[ 

3.3.1. The A phase of^He The A phase of ^He is phenomenologically described by 
Hamiltonian Q, which satisfies particle-hole symmetry C~^'H{—'k)C = — H(k) with 
C = ICtx. Time-reversal symmetry, however, is broken, because the superconducting 
order parameter of Eq. Q is complex. Hence, since = -1-1, Hamiltonian Q 
belongs to symmetry class D. We infer from Table that the Weyl nodes of the 
Hamiltonian Q, which have codimension p = 3 and occur off high symmetry points, 
are protected by a S topological number. In order to derive a formula for this 
topological number it is convenient to rewrite Hamiltonian Q as 

■H(k) = N(k) • T, (13) 

i.e., a dot product between the pseudospin vector N(k) = {hCk), Aok^/kp, Aoky/kp) 
and the vector of Pauli matrices t — (Tx,Ty,Tz). The unit vector field nt = 
N(k)/|N(k)| exhibits singular points at the Weyl nodes K± of H(k). These point 
singularities realize (anti)hedgehog defects in momentum space and are characterized 
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by the Chern number [251 Eg EH EH 
1 


Nc = 


47r 


d k nk • [Sfe^nk x i^fe^nk] 


(14) 


where C denotes a two-dimensional surface parametrized by k = (ki,k 2 ), which 
encloses one of the two Weyl nodes. Choosing C to be a sphere centered at one of 


the two Weyl points, the integral (14) can be straightforwardly evaluated. The nodes 
at K± = (0,0, ±kp) have topological charge Ns^ = ±1. We note that nk restricted to 
the sphere in momentum space defines a map from k S 5^ to the space nk G S^. 


The Chern number Eq. (14), distinguishes between different homotopy classes 
of these maps (BS]. Since the homotopy group 7r2(«S'^) equals the group of the integer 
numbers Z, there is an infinite number of homotopy classes of maps k G 5'^ i—nk and 
Ns 2 can in general take on any integer value. 

Alternatively, instead of considering a sphere 5'^ one can take C to be a plane 
perpendicular to the k^ axis in the Brillouin zone. This yields Nk^ = -1-1 for \kz\ < kp 
and zero otherwise. That is, the unit vector nk for hxed \kz\ < kp possesses a unit 
skyrmion texture in momentum space and measures its skyrmion number. In 


this picture, we can give the topological invariant (14) a more physical meaning by 


reexpressing it in terms of the Berry curvature F(k) = Vk x A(k), i.e., 
1 


Ar = 


27r 


dkxdky Fz(k). 


(15) 


Here, A(k) = i (M_(k)| Vk |w-(k)) denotes the Berry connection (also known as Berry 
vector potential), which is obtained from the negative-energy wavefunctions U-(k) 
of the BdG Hamiltonian. For Hamiltonian ([^, the three components of the Berry 
curvature are given by 


Fx{k) = 


Fyik) = 




F.(k) = 



(U2 1^2 >|2 1 ^0 p (1.2 1 Jc2\ 

3 5 

2 


Alkj-kykz 



(1^2 p \2 1 ^0 p (U2 1 ^2^ 

{h. ^2 ^ '^y) 

3 5 

2 


A2jL2 (1^2 ^2 u2 l2 a 

^x ^y ^F) 



(16) 


2^2 


{k^-klY + ^kl{kl + ki) 


The Berry curvature is plotted in Fig.j^a). We observe from Eq. (16) that the Weyl 
nodes are sources and drains of Berry flux. That is, they act as unit (anti-)monopoles 
of the Berry curvature F(k) and there is a Berry flux of 27r flowing from one Weyl 
node to the other along the z direction. Along the x and y directions, on the other 
hand, the Berry flux is vanishing since Fx(y){'k) are odd in kx (ky) and kz- We also 
note that the Berry curvature is sharply peaked within a shell of width Aq about the 
Fermi surface. 

Due to the bulk-boundary correspondence, the band topology of Weyl 
superconductors gives rise to arc surface states. Consider for example Hamiltonian ([^ 
in a three-dimensional slab geometry with a (100) surface. Because of translation 
invariance along the surface, the three-dimensional (100) slab can be decomposed into 
a family of two-dimensional (10) ribbons parametrized by kz- Ribbons with a kz in 
between the pair of Weyl nodes have a non-zero Chern number, i.e., = -|-1 for 
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Figure 2. Berry curvature in (a) the (px-\-ipy)-w&ve and (b) the {d^ 2 _y 2 -\-idxy)- 
wave Weyl superconductor in the kx-kz plane. The arrows indicate the direction 
of F(k) while the colour scale gives arctan(|F(k)|). We choose Aq = 0.5p; 
realistic values of Aq <C p give a curvature which is much more sharply peaked 
at the Fermi surface. 


\kz\ < kp, while for > kp the topology is trivial. Each ribbon with \kz\ < kp 
can be interpreted as a two-dimensional fully gapped {px + ipy)-weLve superconductor 
with a chiral Majorana edge mode. Hence, the Majorana surface states of the 
Weyl superconductor form a one-dimensional open arc in the surface Brillouin zone, 
connecting the projected Weyl nodes of the bulk gap, see Fig. [^a) and Sec. 4.2 


These chiral surface states give rise to anomalous spin and thermal Hall effects which 


is proportional to the separation of the Weyl nodes in momentum space, see Sec. 5.2 


3.3.2. Chiral (dx^-y^ idxy)-wave superconductor The spin-singlet chiral {dx 2 -y 2 ± 
idxy)-wa,ve superconductor ([^ satisfies both particle-hole symmetry and SU{2) spin- 
rotation symmetry and therefore belongs to symmetry class C in Table (or 
alternatively to class A if one considers only a t/(l) part of the spin-rotation 
symmetry). Hence, as indicated by the classification table, the Weyl nodes of 
Hamiltonian ^ (which occur off high symmetry points) are protected by a 2Z 
invariant, which takes o n on ly even values. This topological number is given by the 
same expression as Eq. (14) but with N(k) = (h(k), Ao(fc^ — ky)/kp,2Aokxky/kp). 
As opposed to the previous example, the unit vector iik = N(k)/|N(k)| for fixed 
\kz\ < kp now describes a double Skyrmion texture with Nk^ = ±2. That is, the 
double Weyl points at K± = (0,0, ±kp) correspond to double (anti-)monopoles of the 
Berry curvature F(k), which is given by |43j 


A,(k) = 


2Alkxkz{kl 


kl) 


(p _ klY 


Fv{^) = 


2AQkykz{k. 


{kl 


kl? 


kl) 




(fc 


2 - 4)2 


A2 


{kl 


kl? 


(17) 
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The Berry curvature is plotted in Fig.j^b). As a result of the topological non-triviality, 
the d-wave Weyl superconductor (§ supports two spin-degenerate, chirally dispersing 
arc surface states, see Fig.|^b). These surface states carry anomalous spin and thermal 
Hall currents. 


3.3.3. Nodal noncentrosymmetric superconductor Nodal NCSs defined by the BdG 
Hamiltonian H(k) in Eq. ([^ satisfy time-reversal symmetry T~^'H[—'k)T = +'H(k), 
with T = ICtq 0 iay, and particle-hole symmetry = — 'H(k), with 

C = ICtx ® ctq. Here, and cTi are the Pauli matrices operating in particle-hole 
and spin space, respectively. Since T^ = —1 and = -1-1, Hamiltonian ([^ belongs 
to symmetry class Dili. According to Table the topological properties of class 
Dili Hamiltonians are characterized by Z 2 topological numbers. While these binary 
invariants can give rise to protected Majorana cone surface states or arc surface states, 
they do not guarantee the stability of the superconducting nodes m- Hence, one 
might conclude that the superconducting nodes of the NCS ([^ are unstable against 
gap opening. In addition to TRS and PHS, however, H(k) in Eq. ([^ also possesses 
the chiral symmetry S~^'H{k)S = —'H(k), with S = —Tx 0 (Jy. Therefore, H(k) can 
be viewed as a member of symmetry class AHI and, as indicated by Table its line 
nodes are protected by a Z invariant, i.e., a winding number. 

The Z 2 and Z invariants of symmetry class Dili and AHI, respectively, can be 
conveniently expressed in terms of the off-diagonal block (/(k) of the spectral projector. 
Eor Hamiltonian ([^, the matrix q(k) takes the following form pTl [22] 


g(k) = 


1 

2Ai,kA2,1 


^k^k l^kl ) O-Q -I- (AjAj^ |lk| 


|lk| 


• (18) 


In Eq. (18) we have introduced the short-hand notation Ak = X + iAt, Hk 


and AJ = Ai^k ± A 2 ,k, where 


£k + iAg, 


Ai.k = + and A 2 ,k = \/ (4)^ + (^k)^ (19) 

are the positive eigenvalues of 'H(k). The unitary matrix ^(k) S U{2) is constrained 
by time-reversal symmetry as iayq^{—'k) = q{'k)iay. 


Winding number The stability of the line nodes of NCS ([^ is ensured by the 
conservation of the winding number [191 HIl IMl EZl ES] 


lAc = ^ j^dfc/Tr[gl(k)afc,g(k)], (20) 

where £ is a closed one-dimensional contour parametrized by fc;, which interlinks 
with a line node. Mathematically speaking, Eq. (20) labels the homotopy classes of 
mappings from £ ~ —)■ q{k) S U{2). Since the first homotopy group of U{2) equals 
the group of the integer numbers Z |66] . i.e., 7ri[f7(2)] = Z, the winding number Wc 
can in principle take on any integer value. Using Eq. we find that Wc for the 

BdG Hamiltonian ([^ can be rewritten in the simpler form 




dki dki [arg (^+ -H iA+) + arg (4 + iA^)] ■ 


( 21 ) 
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Assuming that the superconducting gaps of Hamiltonian ([^ are nonzero only in 
the vicinity of the Fermi surfaces (which is the case for weak-pairing superconductors), 


Eq. (21) can be further simplified to 

1 


lEc = - 


2 ^ 

y—:t koG5' 




( 22 ) 


The second sum in Eq. (22) is over the set of points given by the intersection 
of the helicity-i^ Fermi surface with the one-dimensional contour C. Eq. \22\ 
demonstrates that for a weak-pairing nodal superconductor of the from ([^, the 
topological characteristics of the nodal structure is entirely determined by the phases 
of the gaps in the neighborhood of the positive- and negative-helicity Fermi surfaces. 
This is in complete agreement with the discussion of Sec. |4.3| 

By the bulk-boundary correspondence [53], a non-zero value of Wc signals the 
appearance of zero-energy states at the surface of the NCS. To make this more explicit, 
let us consider a contour C along a line parallel to the (Imn) direction in the bulk 
Brillouin zone, see Fig. With this choice we obtain for Eq. 


W^(Zmn) (k||) = - ^ E (^kj - Sgn ( 

y—± 


AC 


(23) 


where ky denotes the momentum perpendicular to the (Imn) direction, and k± = 
(k||,fc_L^±) and k(j_ = (k||,fc(|_ j_) are momenta on the positive- and negative-helicity 
Fermi surfaces which project onto ky, where k±^± and k'j_ ^ have opposite signs. With 
this, it follows from the bulk-boundary correspondence that zero-energy states appear 
at the {Imn) face of the NCS whenever lF(;mn)(kji) 7 ^ 0. This corresponds to regions 
of the surface Brillouin zone that are bounded by the projected nodal rings of the bulk 
gap, see Fig. [^ In other words, the non-trivial topology of the nodal rings in the bulk 
leads to protected flat-band surface states (cf. Fig. [^c) and Sec. |4.3.3| ). 


Z 2 topological invariant The one-dimensional (two-dimensional) Z 2 invariant of 
symmetry class DIB is given by mEHEollMl ED] 


^ ^ TT Pf [iayq^{K)] 

^ Vdet [icrygT(K)]’ 


(24) 


where Pf is the Pfafhan and the product is over the two (four) time-reversal invariant 
momenta K of the one-dimensional (two-dimensional) centrosymmetric contour C. 
In order for Nc to be well defined, the contour C needs to be centrosymmetric 
(i.e., mapped onto itself under k —>■ —k) and must not intersect with the nodes of 
■H(k). Inspection of Eq. (24) reveals that Nc can only take on the values ±1, where 
— 1 (-1-1) indicates a topologically nontrivial (trivial) character. For weak-pairing 
superconductors the gaps A^ are non-zero only in a small neighborhood of the Fermi 
surface and hence Eq. ^2^ can be simplified to |22| 

= sgn(Aj) sgn(A£) , (25) 


where sgn(A^) represents the sign of the gap at the intersection of C with the 
positive/negative helicity Fermi surface. By the bulk-boundary correspondence |64] . 
Nc = —1 for a one-dimensional contour C indicates the appearance of a helical 
Majorana cone at a surface of the NCS, see Fig.[^a). Similarly, for a two-dimensional 
contour £, Nc = — 1 gives rise to an arc surface state, see Fig. Sb). 
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2D surface 

Brillouin zone 3D Brillouin zone 



Figure 3. The relationship of the bulk gap structure to the surface states 
of a nodal topological superconductor. The left part of the figure shows the 
surface Brillouin zone with the projected Fermi surface indicated in black. Flat- 
band surface states occur within the two regions bounded by the projected nodal 
lines (dark blue and light gray areas). Within these two regions the winding 
number ( |23[ l takes on the values VF = ±1, while outside these regions it is zero; 
the bulk Hamiltonians restricted to a surface momentum in these regions (red 
lines in 3D Brillouin zone) belong to symmetry class AIII and are topologically 
trivial and nontrivial, respectively. The right part of the figure shows the three- 
dimensional bulk Brillouin zone with a spherical Fermi surface (black ellipse) and 
two nodal rings (solid and dashed blue ellipses). 


3.4- Crystalline topological nodal superconductors 

In the previous subsections we have discussed the topological classification of nodal 
superconductors in terms of internal symmetries, i.e., global symmetries that act 
locally in position space (see Table [^. Recently it became apparent that besides 
global symmetries, also crystal symmetries which operate non-locally in position 
space can lead to the protection of superconducting nodes [SH HOI ED ED ED ED- 
Indeed, a complete topological classification of reflection-symmetry-protected nodal 
superconductors was recently obtained in Ref. [60] . Moreover, it has been argued that 
the heavy fermion superconductor UPta can be viewed as a crystalline topological 
nodal superconductor, since its topological surface states are protected by mirror 
symmetry m- 

4. Bound state wavefunctions 

In this section we construct the wavefunctions of the surface states of a topological 
superconductor within an approximate low-energy effective theory. Such so-called 
“quasiclassical” methods are a standard technique in the study of inhomogeneous 
superconductivity ED ED Ei ED- The focus here is to demonstrate how this approach 
can be used to explore the surface physics of topological superconductors. 
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^.1. Deriving the low-energy theory 


We consider a surface of a topological superconductor, with normal defining the axis 
r±. The superconductor is assumed to occupy the half-space r± > 0, and translational 
invariance parallel to the surface ensures that the transverse momentum ky is a good 
quantum number. Our aim is to solve the BdG equation 

#^(ki|;ri) =i;^'(k,|;ri) (26) 

for states with energy comparable to the superconducting gap. The Hamiltonian 
i?(kj|) appearing here is related to the bulk Hamiltonian H(k) by 




(27) 


The BdG equations can be solved numerically, e.g. by exact diagonalization of a lattice 
regularization of the Hamiltonian. Superconductivity is a low-energy phenomenon, 
however, where the characteristic energy scale of the gap A is typically much smaller 
than the Fermi energy Ep, and it is at such low energy scales that the surface states 


exist. Not only is it computationally advantageous to restrict our solution of Eq. (26) 


to the low-energy sector, but it also permits at least semi-analytic solutions which 
provide a more transparent understanding of the physics. This is the motivation for 
the quasiclassical theory of superconductivity. A full development of this theory, 
in particular the construction of the quasiclassical Green’s function and the self- 
consistent determination of the pairing potential, is rather involved and beyond the 
scope of this review. We therefore sketch the theory below, and refer interested readers 
to the extensive literature on the subject, e.g. [761 [771 [781 [79]. 

The first step in constructing the effective theory is to bring the noninteracting 
Hamiltonian /i(k) into a diagonal form. The low-energy quasiparticle excitations in 
the superconductor occur within a narrow momentum shell about the normal-state 
Fermi surface of width ^ A/S|vf| ^ |ki?| where vp is the Fermi velocity. The 
Fermi velocity is regarded as constant across this momentum shell, which allows us 
to linearize the low-energy Hamiltonian [/i(k)]a^Q, « hxp^a • (k — ki^’^a) where kj’^a lies 
on the Fermi surface of the a band. We now examine the pairing matrix A(k). The 
superconducting gap is taken to be constant along the direction of the Fermi velocity at 
each point on the Fermi surface, i.e. A(k) « A(kF^Q,). The solution of the low-energy 
effective theory is therefore only sensitive to the value of the superconducting gap in 
the immediate vicinity of the Fermi surface. Pairing between states on nondegenerate 
Fermi surfaces is not considered. 

An appropriate ansatz for an eigenstate of the low-energy Hamiltonian is 


^(k|i;r) =51' 


u. 


phk||T||-|-[fe3+5fej]r-x) 


(28) 


where the sum j is over all scattering channels, Oj are C-number coefficients, the 
vector kj = (ky, kj) lies on the Fermi surface, \Skj \ ^ \kj \ is a small energy-dependent 
deviation from the Fermi wavevector k^, and the multi-component BCS coherence 


factors Uj and Vj 


are functions only of kj 


and the energy. The coefficients aj 
in Eq. (281 are determined by imposing boundary conditions on the wavefunction 


at the surface, which can in principle be derived from the underlying microscopic 
theory. This presents the greatest challenge faced by the low-energy effective theory, 
as there is no general prescription of how this should be done, especially for multiband 
systems or those with strong spin-orbital mixing. 
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In order to avoid this complication, in the following we restrict our attention to 
a commonly-used microscopic model of the normal state where boundary conditions 

can be explicitly and easily derived (TJ [TU [T^ [521113113111331 @311331 EZl [7^ ISHIM1 @31 
IMl 1551155] . Specifically, we adopt the noninteracting Hamiltonian 




do + Aik • (T -|- Vwork0(?'_L) 


(29) 


This describes a free electron gas, including antisymmetric spin-orbit coupling of 
strength A so as to model both centrosymmetric (A = 0) and noncentrosymmetric 
(A ^ 0) systems. In the latter case, the spin-orbit vector Ik has the following 
continuum-limit representations for the different point groups 

{aikx + a2ky)yi N {a^k^ + aiky)y + a^k^z, 
kyy. - k^y, 


Ik 


kx{i+g2[kl + 


kyil + g 2 [kl + kl])y + fc^(I -f g 2 [kl + kl])z, 


kx{kl - kl)yi + kyikl - kl)y + k^{kl - kl )2 


C2 

o’'(30) 

Td 

The coefficients and g 2 appearing in the expressions for the C 2 and O point 

groups are arbitrary. Our Hamiltonian also includes a confining potential Kvork for 
the surface which we identify with the work function. For our low-energy theory it is 
reasonable to take the limit Kvork —>■ 00 , and so we implicitly include this potential as 
the boundary condition that the wavefunctions vanish at the surface. 


4-2. Bound states of Weyl superconductors 

As an illustration of the quasiclassical method, we consider the subgap bound 
states appearing at the (100) surface of the Weyl superconductor models introduced 
in Sec. |2.Il Take a point ky = {ky,kz) in the two-dimensional surface Brillouin 
zone that lies within the projection of the Fermi surface, i.e. |kj|| < kp. This 
surface wavevector corresponds to two points on the Fermi surface p = (fc,j,,kjj) and 

p' = (— k||), where k^ = ^Jkf.— |kj|p. The bound state wavefunction is then a 
superposition of evanescent states at these momenta 


4'(k|i;r) = a 


^-(ik^+K.-p,)x I ik 


7p 7 V 7p' 

where Ok are C-number coefficients, and the coherence factors 7k and inverse decay 
length Kk are defined 

I 

Ak L 

m^|Akp - E'^ 
h?\kx\ 


7k = 


Kk = 


isgn(/c3;)y/|Akp - A2 


(31) 


(32) 


(33) 


The ansatz Eq. (311 satisfies the criteria that the wavefunction vanishes in the bulk. 


The requirement that the wavefunction also vanishes at the surface leads to the 
condition on the coherence factors 

7p - 7p' = 0 . (34) 

This can be analytically solved for the surface bound state energies of the p-wave and 
d-wave Weyl superconductor: 

/ - \ 1/2 _ 

Ao (1 - kf\ ky px; + ipy 


Em) = 


(i-t?)( 


1 - |k||P - kl]sgn{ky) 


^dxy 


(35) 
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Figure 4. Topological surface states in Weyl superconductors. The surface 
bound-state spectra is plotted as a function of the (100) surface momentum 
{ky,kz) for (a) the (px + 2 py)-wave and (b) the {d^ 2 _y 2 + idxy)-'wsive Weyl 
superconductor. The energy of the bound states is indicated by the color scale, 
where the jump from cyan to magenta occurs at zero energy; light grey indicates 
the absence of a surface state. The light green line shows the projection of the 
Fermi surface. 


where = k^/kp- 


The bound state spectrum is plotted in Fig. In both cases, 
we observe zero-energy arc surface states running between the nodes at the north 
and south poles of the Fermi surface. Whereas there is only a single such arc state 
in the p-wave Weyl superconductor, there are two arc states in the d-wave Weyl 
superconductor, consistent with the Chern numbers Eq. (15) of Nk^ = 1 and Nk^ = 2 
for \kz\ < kp, respectively. 


4-3. Bound states of NCS 


We now turn our attention to the surface bound states of a NCS. We restrict 
ourselves to the limit Ep ^ Amax{|lk|} ^ |A|, so that we can ignore the spin¬ 
splitting of the Fermi surfaces in the low-energy theory, but the triplet d-vector is 
still constrained to be parallel to the spin-orbit coupling vector, i.e. dk = AJk, 
and hence the physical description outlined in Sec. |2.2| is still applicable. This is 
a common approximation [53J [731 IHSl ISZl |HH1 IMl [301 IHI] as accounting for the 
spin-orbit splitting of the Fermi surfaces complicates the analysis but generally only 
quantitatively modifies the surface state spectra. 

A similar (but more involved) argument to that for the Weyl superconductors 
above gives the implicit equation [331131173] for the bound states energies 

0 = (7p' -7p)(7p' -7p)(|lpl|lp'l - Ip • Ip') 

+ (7p/ -7p)(7p' -7p)(|lpl|lp'l +lp • Ip')- (36) 

where 



E - i sgn(wF.u (k))y|A^|2 - 


(37) 


where Ui7',_L(k) is the component of the Fermi velocity in the direction normal 
to the boundary. Bound-state energies must lie within the minimum gap, i.e. 
|E| < min{|Ap I, |Ap,|}. Although the general solution must be obtained numerically, 
analytic conditions can be derived for the existence of zero-energy states. Substituting 
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Figure 5. Helical Majorana and arc states in NCS. The surface bound-state 
spectra is plotted as a function of surface momentum ky for: (a) (110) face of C 2 
NCS with r = 0.4 and = 1; (b) (100) face of NCS with r = 0.25; 

(c) (100) face of O NCS with r = 0.1 and g 2 = —1.5; (d) (110) face of NCS 
with r = 0.2. The energy of the bound states is indicated by the color scale, 
where black corresponds to zero energy and yellow to the maximum energy F^max; 
light grey indicates the absence of a surface state. The light green line shows 
the projection of the Fermi surface. As discussed in the text, in each case the 
zero-energy state at ky = 0 is a topological helical Majorana. The C 2 and Cav 
NCS also display topologically-protected arc states. 


the zero-energy form 


of 7k 


= —i sgn(t)i?_j_(k)A^) into Eq. (36) we obtain five 


conditions for zero-energy states in terms of the spin-orbit coupling vector and the 
gap signs [22] : 


(i) The spin-orbit vectors Ip and Ip' are parallel, and the gap signs satisfy sgn(A+) = 
sgn(Ap,) sgn(A+) = sgn(Ap). 

(ii) The spin-orbit vectors Ip and Ip' are antiparallel, and the gap signs satisfy 
sgn(A+) = sgn(A+) ^ sgn(Ap,) = sgn(Ap). 

(iii) The gap signs satisfy sgn(A+) = sgn(Ap,) and sgn(Ap,) ^ sgn(Ap). 

(iv) The gap signs satisfy sgn(A+) ^ sgn(Ap,) and sgn(Ap,) = sgn(Ap). 

(v) The gap signs satisfy sgn(A+) ^ sgn(A+), sgn(Ap,) ^ sgn(Ap), sgn(Ap,) ^ 
sgn(A+), and sgn(A+) sgn(Ap). 


These Hve conditions for a zero-energy state include the topological criteria defined 
in Sec. 3.3.3| Only a state satisfying the first condition is not topological: It cannot 
be a Kramers-degenerate Majorana mode as the antisymmetric spin-orbit vectors at 
p and p' are not parallel if ky lies at a time-reversal invariant point in the surface 
Brillouin zone; Nor can it be an arc surface state protected by a Z 2 number, since the 
condition on the gap signs requires that there be a node in the plane defined by p 
and p'. On the other hand, since condition (ii) requires different sign of the positive 
and negative helicity gaps, it is clearly consistent with the definition of the topological 
number Eq. (25) for both these topological states. It also describes nontopological 
states, however, as these conditions are not sufficiently restrictive. The conditions 
(iii)-(v) exclusively describe flat band states with nontrivial winding number, and 
they are equivalent to the topological criteria Eq. (23). Specifically, conditions (iii) 
and (iv), which correspond to a situation where only one of the helical gaps has a sign 
change between the two sides of the Fermi surface, describe topologically protected 
nondegenerate zero-energy states with winding number ITpm„) = ±1. In contrast. 
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Figure 6. Topological fiat-band states in nodal NCS. In the top row we plot 
the surface bound-state spectrum as a function of surface momentum ky for (a) 
the (100) surface of a C 2 NCS with r = 0.25 and = 1, (b) the (101) 

face of a C 4 V NCS with r = 0.25, (c) the (111) face of O NCS with r = 0.25 and 
g 2 = —1.5, (d) and the (111) face of NCS with r = 0.15. The colour scale is 
the same as in Fig. In the bottom row we plot the corresponding variation of 
the VF(Zmn)(^11 )• Dark blue (gray) corresponds to (^||) = +1 ("l)^ while 

light blue represents. VF(^^^)(k||) = 0. Figure (a), (b), (c) adapted from |83| . 


(v) requires sign changes of both helical gaps and gives doubly degenerate states of 
winding number W(^imn) = i2. 

4-3.1. Helical Majorana Helical Majorana states appearing at the centre of the 
surface Brillouin zone (ky = 0) are relatively common for majority-triplet NCS at 
high-symmetry surfaces. In each of the examples shown in Fig. the zero-energy 
state at the zone centre has this topological protection. The other zero-energy states 
are present elsewhere in the surface Brillouin zone in each of these examples, on the 
other hand, do not satisfy the necessary criterion. 

4 . 3 . 2 . Arc states Topologically protected arc states are rare, as the condition that 
Ik be symmetric upon reflection about some plane places severe constraints on the 
point group symmetry. There are two example in Fig. portions of the arc states 
at the (110) surface of the C 2 V NCS shown in panel (a), and the entire line of zero 
energy states connecting the projection of the nodal rings at the (100) surface of the 
C 4 V NCS. In contrast, the arc states at the (110) surface of the Td NCS [panel (d)] 
do not have topological protection, as the topological condition only holds for planes 
which intersect gap nodes. 

4 . 3 . 3 . Flat-band states Flat-band states are the most generic topological surface 
state of nodal NCSs, appearing at any surface where the projection of the nodal rings 
encloses a finite area and do not exactly lie on top of another. We illustrate this with 
an example from each point group in Fig.[^ In the top row we plot the spectrum in the 
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Figure 7. Spin polarization of the surface states at the (101) surface of a Cav 
point group NCS with r = 0.25. Panels (a) and (b) show the spin polarization 
along the x- and y- axes of the crystal, respectively. The 2 :-spin polarization 
is vanishing. The energy of the states is shown in Fig. ib). Figure adapted 
from [83]. 


surface Brillouin zone, while in the bottom row we show the corresponding variation of 
the winding number lV(/mn)(k||)- The regions where the topological number evaluates 
to ±1 are bounded by the projected line nodes of the gap, and exactly matches the 
location of the zero-energy flat bands in the surface spectrum. By the bulk-boundary 
correspondence this ensures the topological protection of the zero-energy flat bands 
and implies that these states are nondegenerate. 


4-3.4- Spin polarisation of surface states The surface bound states of NCSs typically 
display a strong spin polarization, which is a consequence of the broken inversion 
symmetry. The /r-component of the spin polarization of a bound state at transverse 
momentum ky and energy E is defined |83j as the expectation value 

p'^(ii;,kj|)=^“drxd/t(k„;r^)(^ ^ vl/(k||; rj.) (38) 


where 'I'(kjj;r_L) is the wavefunction. Due to the symmetries of the Hamiltonian 
(time-reversal, particle-hole, and chiral), the polarization of the edge states is even 
in energy and odd in the surface momentum, i.e. p^(i?,ky) = pf*)—i?,k||) = 
—p^{E, —ky), and there is no spin accumulation at the surface. The spin polarization 
is sometimes defined only with reference to the electron-like components of the 
wavefunction [E1IHH1I57], but the definition Eq. is more reasonable as this quantity 
couples to an applied Zeeman field. 

The spin polarization of the edge states at the (101) surfaces of a C^y NCS is shown 
in Fig. Both the dispersing and the flat-band edge states display a pronounced spin 
texture. The spin polarization does not follow from the topology, however, but rather 
arises from the unequal contribution of positive- and negative-helicity states to the 
wavefunction. The spin polarization of the surface states must match that of the 
continuum states where the two intersect, in particular at the bounding nodes of 
the flat band states. Despite its nontopological origin, the spin polarization of the 
topological flat-band states is exploited in several proposed experimental tests of their 
existence. 
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5. Experimental signatures of topological surface states 

Directly detecting the surface states of a topological superconductor requires 
measurement techniques sensitive to the surface properties. This is a much greater 
experimental challenge than demonstrating the existence of gap nodes in the bulk, 
as it generally requires the preparation of high-quality surfaces or the incorporation 
of the topological superconductor into a heterostructure device. These hurdles are 
not insurmountable, however, as evidenced by the observation of surface states in 
unconventional superconductors such as the cuprates [HlllSSl lU and Sr 2 Ru 04 [33]. 
In this section we examine different experimental signatures of topological surface 
states. 

5.1. Tunneling conductance 

Conductance spectroscopy experiments represent the dominant approach to probing 
the surface states of unconventional superconductors [771 Hg. Since the conductance 
is sensitive to the surface density of states of the superconductor, it can provide 
direct evidence of surface states. A number of authors have studied the conductance 
signatures of topologically protected surface states of NCS 1211123 IMm 113 Ha IMl 

IWIIMIIMI [Tool 188113T]. 

The simplest model of a conductance experiment consists of a junction between 
a topological superconductor and a normal metal with an insulating barrier at their 
interface. The conductance as{E) at zero temperature is given by a generalization of 
the familiar Blonder-Tinkham-Klapwijk formula [IM llOlj . 

(39) 

where and b'^’^ are the Andreev and normal reflection coefficients, respectively, 
for spin-cr electrons injected into the superconductor at interface momentum ky. In the 
presence of spin-orbit coupling, both spin-preserving (tr' = a) and spin-flip {a' = —a) 
Andreev and normal reflection coefficients are included. The scattering coefficients 
are determined by solution of the BdG equation for an electron at bias energy E = eV 
injected into the superconductor from the normal lead. The insulating tunneling 
barrier is commonly modeled as a (5-function with strength U = ZEp/kp [77] ■ The 
value of the dimensionless barrier parameter Z defines different experimental regimes: 
small values of Z < 1 are least sensitive to surface states, and are typically used to 
model point contact Andreev spectroscopy experiments m; larger values Z >1 give 
clear signatures of the surface states and are characteristic of tunnel junctions and 
scanning tunneling microscopy. We will concentrate upon the latter regime. 

Tunneling conductance spectra at the (100) and (101) interfaces of a C 4 „ NCS are 
shown in Fig. While the majority-singlet cases (r > 0.5) show a clear gap feature, 
majority-triplet pairing (r < 0.5) is characterized by substantial subgap conductance 
due to resonant tunneling through the surface states. At the (100) interface, the 
dispersing arc states give a broad feature centred at zero bias which completely fills 
the gap. This is also present at the (101) interface, but here the most striking aspect is 
the sharp zero-bias conductance peak due to tunneling into the zero-energy flat-band 
states [77]. The sharp zero-bias peak is a robust indicator of the flat-bands, as it 
arises from their singular contribution to the surface density of states. On the other 
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Figure 8. Tunneling-conductance spectra at the (100) and (101) interfaces of a 
C 4 ,v NCS with tunneling barrier strength Z = 4. The results are normalized by 
the normal state conductance in the absence of superconductivity. In both 
panels the singlet-triplet parameter r is tuned from r = 0 (purely triplet) at the 
top to r = 1 (purely singlet) at the bottom in steps of 0.05. The curves at r = 0, 
0.25, 0.5, 0.75, and 1 are highlighted in black. The conductance spectra are offset 
by 0.2 for clarity. 


hand, fine structures in the dispersing states gives rise to a much richer variation in 
their contribution to the conductance, and there is no universal signature [13IZ3- In 
particular, there is no way to distinguish the dispersing surface states of a NCS from 
those of a Weyl superconductor using simple tunneling conductance measurements. 

Many variations on the basic junction geometry are possible, e.g., replacing 
the normal lead by a metallic ferromagnet pTl fTODllTUH) . or assuming a spin-active 
barrier |83j . This allows conductance measurements to probe the spin structure and, 
indirectly, the degeneracy of the surface states. Applying a Zeeman field to the sample 
is predicted to shift the spin-polarized flat bands away from zero energy, thus splitting 
the zero-bias peak PH 1104) . The absence of splitting for certain field orientations 
would strongly suggest nondegenerate bands, as doubly-degenerate bands are expected 
to be spin degenerate and split for any orientation of the applied field |105) . On the 
other hand, the stability of nodal NCSs in such Zeeman fields is uncertain |l()b[ 1107] . 
A less invasive variation on this approach is to replace the insulating tunneling barrier 
with a ferromagnetic insulator |1051 [83] . 

The topological classification of nodal superconductors is only approximately valid 
in the presence of disorder. Even for a pristine bulk system, however, disorder at the 
surface or at the interface in a tunneling junction is inevitable. So long as the disorder 
is not too strong the conductance features discussed above remain detectable, although 
they are broadened by the finite lifetime of the quasiparticles [108111091 IST] 111011111) . 
Extremely disordered surfaces produce “dead” insulating layers, which redefine the 
edee of a suoerconductor |110l 111111112). This is illustrated in Fies. ia) and [^b) , 
which show the layer-resolved density of states of a {dxy + p)-wave NCS [H] on a 
square lattice in the presence and absence of strong edge disorder, respectively. For 
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Figure 9. Layer-resolved local density of states pi of a {dxy -|-p)-wave NCS in a 
ribbon geometry with (01) edges, which displays zero-energy fiat bands. Panel (a) 
corresponds to a dirty NCS with non-magnetic Gaussian disorder in the outermost 
layer of strength 7imp = 200At, while panel (b) represents a clean NCS. In the 
former, the outermost layer is “dead” and the topologically protected edge states 
develop in the second layer. Figure adapted from Ref. [no]. 


very strong edge impurities, the states in the outermost layer are almost fully localized, 
while in the second and third inward layers new weakly disordered states reemerge. 

5.2. Surface currents 

A key signature of the chiral surface states of Weyl superconductors is that 
they show anomalous spin and thermal Hall conductance in the direction of their 
propagation [JHl 1113] . Although the zero-temperature Hall conductances are directly 
related to the Chern number, they also depend upon the number of chiral edge modes 
in the system, and so they are proportional to the separation of the Weyl points 
in momentum space. This is in contrast to the universal quantization in a fully- 
gapped time-reversal-symmetry breaking topological superconductor |114l I115L1116] . 
Furthermore, the spin Hall conductance in the triplet Weyl superconductor depends 
on the orientation of the magnetic field with respect to the d-vector, as not all spin 
directions are conserved m- While the measurement of the anomalous spin Hall 
conductance is likely to be very difficult |114j , the anomalous thermal Hall conductance 
is predicted to be experimentally accessible for realistic materials [13J 1117] . 

The chiral surface states may also carry a surface charge current, but it is 
not topologically protected since charge is not conserved in a superconductor. In 
particular, its detection is extremely difficult as it will be compensated by Meissner 
screening currents |118j . While the current carried by the bound states is localized 
within a few coherence lengths ^ of the surface, however, the screening currents build 
up on the scale of the penetration depth A, and so in extreme type-H superconductors 
a finite surface current density is expected. Precision magnetometry has so far failed to 
detect surface charge currents in Sr 2 Ru 04 |119) . however, despite otherwise compelling 
evidence that this is a chiral p-wave superconductor |120L I121j . The charge current is 
not universally related to the Chern number, however, and even ignoring the Meissner 
effect, it depends sensitively upon microscopic details of the superconductor |122[I123] . 

Surface currents might also arise in time-reversal invariant NCS. Ignoring spin- 
orbit coupling, the simplest model of a majority-triplet C 4 v NCS is a superconducting 
analogue of the quantum spin Hall insulator |20] , and the helical edge states therefore 
carry a surface spin current [HH] II241 151] . On the other hand, a surface spin current 
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Figure 10. Fourier-transformed quasiparticle interference patterns for arc states 
at the (010) surface of a point group NCS. Results for (a) nonmagnetic 
scatterers and (b) magnetic scatterers with impurity moment along the a:-axis. 
The weak and diffuse features in (a) reflect the suppression of backscattering due 
to the spin-polarization of the edge states; backscattering occurs for magnetic 
scattering, giving rise to the much stronger and sharper features in (b). Figure 
adapted from Ref. [65]. 


carried by states above the gap is present in the nontopological phase [57], and the 
spin current is strongly suppressed by the spin-orbit coupling [55] . Charge currents are 
also predicted to arise at the interface between NCSs and ferromagnets [5T1 11251 [l2hj . 
In particular, at surfaces with spin-polarized flat bands, the coupling to the exchange 
field of the ferromagnet generates a chiral electronic structure at the interface. The 
resulting charge current is very pronounced at low temperatures and shows nonanalytic 
dependence on the exchange field strength [51] I125j . 


5.3. Other proposals 


There are a number of proposals for evidencing chiral or helical edge states based upon 
interference effects in the magneto conductance of superconductor loops [12711128| , or 
nonlocal conductance at domain walls in Weyl superconductors |129| . Although such 
experiments could give unambiguous demonstration of the edge states, these proposals 
assume a fully-gapped pairing state, so it is unclear if they would apply to a nodal 
system. 

The surface states of a topological superconductor could be more directly 
demonstrated using quasiparticle interference spectroscopy [55lll30| . which has already 
been successfully used to evidence topological insulator surface states m- This 
can probe the spin character of the surface states, as scattering between states with 
opposite spin polarization (i.e. backscattering processes) are suppressed, unless the 
scattering is by magnetic impurities. Predictions [65] for NCS arc states are shown in 
Fig. 10 whereas nonmagnetic scattering gives a weak and diffuse signal [panel (a)], 
magnetic scattering produces a strong and well-defined image of the arc states [panel 

(b)]. 


6. Survey of candidate materials for nodal topological superconductivity 

In certain materials with strong electronic correlations, topological superconductivity 
can occur naturally due to unconventional m, be. non-phonon-mediated, pairing 
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interactions. These unconventional pairing mechanisms (e.g., due to spin fluctuations) 
are in general quite anisotropic and hence give rise to nodal structures of the 
superconducting gap. As a consequence, topological superconductors with nodes are 
far more common than those with a full excitation gap. An important exception to this 
rule is the fully gapped spin-triplet superconductor Sr 2 Ru 04 |12()[I121) . In this section 
we give a brief overview of candidate materials for nodal topological superconductivity 
with either spin-triplet or spin-singlet pairing, which we list in Table and Table 
respectively. We note that due to strong atomic spin-orbit coupling in many of these 
materials, the spin symmetry of the Cooper pairs is determined by the pseudospins 
of quasiparticles with opposite momenta, rather than by the real spins of electrons. 
Thus, the pairing symmetries may be better distinguished in terms of the parity of the 
orbital pair wave function; for simplicity, we nevertheless refer to odd- and even-parity 
states as spin triplet and singlet, respectively. 


6.1. Nodal superconductors with spin-triplet pairing 


A nodal spin-triplet pairing state has been established unambiguously in superfluid 
^He [53l El]. That is, the A phase of ^He is known to be described by the Anderson- 
Brinkman-Morel (ABM) state [52], a three-dimensional superfluid with nodal points 
at the north and south poles of the Fermi sphere. The ABM state is a member 
of symmetry class D with codimension p = 3 (cf. Sec. 3.3.1). Hence, according to 
the classification of Table the point nodes of ^He-A are topologically stable and 
protected by a two-dimensional Chern number [23 [23 [H] [28] . 

Here, however we want to focus on nodal topological superconductors rather 
than superfluids. A very promising group of materials for nodal topological 
superconductivity are the noncentrosymmetric superconductors (NCSs) [1331 1134] . 
These systems have extremely diverse properties, however, as they are defined by 
their crystal structure as opposed to chemical composition. In particular, the parity¬ 
mixing of the gap depends upon microscopic details. The most interesting case for us 
is when the spin-triplet and spin-singlet components are comparable in magnitude, 
which generically leads to a nodal superconducting state as discussed in Sec. 
Sizable spin-triplet pairing components can be expected to occur in strongly correlated 
NCSs where the onsite Coulomb repulsion suppresses the spin-singlet s-wave channel, 
and where there are strong noncentrosymmetric spin fluctuations [1351 11361 1137] . 
These conditions are believed to hold in the heavy fermion C 4 y point group NCS 
CePtaSi [138] . leading to an (s -|- p)-wave state [136] . Indeed, there is significant 
evidence of line nodes and an unconventional spin pairing state in this material 
11401114111142] . A similar pairing state seems to be realized in CeIrSia [14311144] and 
CeRhSia [145] . which are also heavy fermion NCS with C 4 „ point group. Not all NCS 


with sizable spin-triplet gap are so strongly correlated: Replacing Pd by Pt in the O 
point group NCS Li 2 (Pdi_ 3 ;Pta;)B appears to tune the system from a conventional s- 
wave superconductor to a nodal triplet superconductor [1461114711148] , but there is no 
evidence that the correlation strength is dramatically altered by this substitution. The 
gap structure of LaNiC 2 {C 2 v point group) is controversial [1491115011151] . Muon spin 
rotation experiments evidence time-reversal-symmetry breaking in LaNiC 2 [149] and 
the related (centrosymmetric) LaNiGa 2 [152] . suggesting nonunitary triplet pairing. 

A number of Uranium-based heavy fermion compounds show strong evidence 
of unconventional triplet pairing states with nontrivial topological properties mB 
EMI EHH EM] . Among them, UPta is probably the most promising candidate for 
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Table 2. List of candidate materials for nodal topological superconductivity 
with (majority) spin-triplet pairing. Note that in some materials the evidence 
is contradictory. NCS: Noncentrosymmetric superconductor. HF: Heavy fermion 
superconductor. FM: Ferromagnetic superconductor. NMR: Nuclear magnetic 
resonance. SH: Specific heat. UA: Ultrasound attenuation. TC: Thermal 
conductivity. PD: London penetration depth. Hc 2 '- Upper critical field 


Material 

Type 

Evidence for 
triplet pairing 

Evidence for nodes 

Probable pair¬ 
ing symmetry 

A phase of ^He 

superfluid 

NMR, magnetiz. 

SH 

chiral 

CePtaSi 

NGS, HF 

indirect 

Pn. NMR. etc. |1 391114011141 11142) 

(s -I- p)-wave 

GelrSigt 

NGS, HF 

NMR [mi \m\ 

NMR 114:j||144] 

(s + p)-wave 

GeRhSigt 

NGS, HF 

h, 2 [m 


? 

Li2Pt3B 

NCS 

NMR [147J 

PD, NMR, SH )T46i niTi rri^ 

(s -I- p)-wave 

LaNiCa 

NCS 

indirect jl49) 

PI) jb^Uj 

nonunit ary 

LaNiGaa 

centro. 

indirect [152] 


nonunit ary 

URhGe 

FM, HF 

indirect |15.3) 

SH [154] 

p-wave 

UGoGe 

FM, HF 

NMR [T55] 

indirect |156) 

p-wave 

UGeat 

FM, HF 

H,2 1157111581 

NMR [159) 

p-wave 

UPta 

HF 

NMR )16()| 

SH, UA, TC 11611 

chiral /-wave 

UBeia 

HF 

NMR )162| 

SH, JNMR 11 63111641 

nodal 


t 


superconducting under pressure 


spin-triplet superconductivity, and has received much attention CM]- This material 
has three different nodal superconducting phases, and there is strong evidence that 
one of these realizes a time-reversal-symmetry breaking chiral /-wave state, which 
is predicted to show Weyl nodes [11711169) . The superconductivity of UBeia is less 
well-understood, but evidence of surface states has been reported |170j . Intriguingly, 
some of these systems, in particular URhGe, UCoGe, and UGe 2 , show an apparently 
cooperative coexistence of ferromagnetism and superconductivity |1651 1168) . That 
is, close to a ferromagnetic quantum critical point ferromagnetic spin fluctuations 
mediate spin-triplet pairing. The natural pairing state for such systems would involve 
equal-spin-pairing Iff) and j//) on the majority- and minority-spin Fermi surfaces, 
respectively. The crystal space group of most ferromagnetic superconductors is of low 
symmetry, which results in strong anisotropies of the magnetic fluctuations, favourable 
to the observed nodal pairing states. Topological properties of these systems are 
largely uninvestigated, although the existence of arc states connected with Weyl nodes 
has been proposed |130) . 

The relative rarity of intrinsic spin-triplet superconductors has lead to a 
number of proposals to engineer such states in heterostructures involving spin-singlet 
superconductors |48l |49l [801 11711 1172) . A key to these proposals is the presence 
of strong spin-orbit coupling, which, under appropriate conditions, can accomplish 
a change in parity of the Gooper pairs, resulting in a sizable proximity-induced 
spin-triplet gap. For example, superlattices of topological insulators and s-wave 
superconductors are proposed to realize a p-wave Weyl superconductor [151 [15] . 
Incorporating high-temperature superconductors into such devices has the great 
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Table 3. List of candidate materials for nodal topological superconductivity 
with (majority) spin-singlet pairing. HF: Heavy fermion superconductor. 
SL: Superlattice. ARPES: Angle-resolved photoemission spectroscopy. STM: 
Scanning tunneling microscopy. NMR: Nuclear magnetic resonance. PD: London 
penetration depth. SH: Specific heat. TC: Thermal conductivity. MT: Magnetic 
torque. PKE: Polar Kerr effect. ^SR: Muon spin rotation. 


Material 

Type 

Time-reversal 

symmetry 

Evidence for nodes 

Probable pair¬ 
ing symmetry 

YBa2Cu306_|_a;, 
La 2 -a;Sra;Cu 04 , etc. 

high-temp. 

supercond. 

Yes 

ARPES, STM, 

NMR, PD, etc. [51 [33 

da,2_j,2-wave 

CeCu2Si2 

HF 

Yes 

indirect |174| 

SH, TC, NMR, 

d-wave 

CeCoIns 

HF 

Yes 

STM I3H1137113311331117hl 

da;2_y2-wave 

Celrins 

HF 

Yes 

SH. TC. NMR [3711173 

d-wave 

CeRhIns't' 

HF 

Yes 

SH |17fi| 

d-wave 

URu2Si2 

HF 

No (MT HU, 
PKE [T77[1 

SH, NMR, TC [13 SI] 

(d ± id)-wave 

SrPtAs 

pnictide 

No (/rSR |46|1 

indirect [^ 

(d ± id)-wave 

CeCoIns/YbCoIns 

SL 


indirect [178] 

likely d-wave 

Cua, (PbSe)5 (Bi2Se3)6 

SL 


SH |179| 

line node 


t 


superconducting under pressure 


advantage that the engineered spin-triplet superconductor could be realized at much 
higher temperatures than their intrinsic counterparts [501 11711 1172) . Particularly 
interesting is proximity-induced pairing in the surface states of a topological 
insulator [13 HU- Recently, cuprate superconductor-topological insulator devices 
have been fabricated m- However, the proximity-induced gap in these devices is 
found to be nodeless. 

6.2. Nodal superconductors with spin-singlet pairing 

High-temperature cuprate superconductors are probably the most extensively studied 
unconventional superconductors |ld2j . While there is no common consensus on the 
pairing mechanism, it is by now generally accepted that the pairing symmetry in the 
cuprates is of dx 2 _y 2 -wave form [53 ES]- That is, the superconducting gap changes 
sign along the Fermi surface, giving rise to four topologically protected point nodes. 
The superconducting state of the cuprates belongs to symmetry class CH in Table 
since it preserves both time-reversal and SU{2) spin-rotation symmetry. Hence, it 
follows from the classification of Table that the point nodes are protected by a 2Z 
topological number, i.e., by a winding number which takes on the values ±2. Due 
to the bulk-boundary correspondence, the topological characteristics of these point 
nodes lead to the appearance of doubly degenerate (i.e., spin-degenerate) flat band 
edge states [33 133 ■ These zero-energy edge modes have been observed in tunneling 
experiments on the cuprate material YBa 2 Cu 306 -a; [92l [93l l94] - 

A superconducting gap with d-wave symmetry also occurs in the Cerium-based 



























CONTENTS 


28 


heavy fermion compounds CeCu 2 Si 2 |180) and CeXIns, where X is a transition 
metal (Co, Ir, Rh) [Ml 133 ■ In all these materials superconductivity appears close 
to an antiferromagnetic quantum critical point and is believed to be mediated by 
spin fluctuations [1671116811132] . In contrast to the ferromagnetic superconductors, 
however, magnetic order in the Cerium compounds competes with superconductivity. 
NMR measurements indicate spin-singlet pairing and in combination with specific 
heat and thermal conductivity measurements yield evidence for line nodes in the 
superconducting gap. Therefore, it seems likely that the order parameter symmetry 
in these heavy fermion superconductors is of d-wave form. This has been confirmed 
for CeCoIns by a careful analysis of the quasiparticle interference patterns observed 
in scanning tunneling spectroscopy [38l [39] . 

The table lists two compounds which show signs of unconventional spin-singlet 
superconductivity with broken time-reversal symmetry. The first one is URu 2 Si 2 m, 
in which superconductivity emerges out of a somewhat mysterious “hidden order 
phase”. Recently, a polar Kerr effect has been observed, which appears with the onset 
of superconductivity Together with specific heat mi and thermal conductivity 

measurements [40], this suggests that URu 2 Si 2 is a chiral d-wave superconductor, 
which breaks time-reversal symmetry. Such a state can support Weyl nodes |43j . 
The second candidate for time-reversal symmetry breaking superconductivity is the 
pnictide material SrPtAs. Muon spin-rotation measurements on SrPtAs have revealed 
time-reversal-symmetry breaking in the superconducting state [46] . This material is 
also predicted to realize a chiral d-wave state with Weyl nodes [33. 

Finally, we comment on artificial superconducting structures. Recently, superlat¬ 
tices of the ( 1^2 _y 2 -wave superconductor CeCoIns and the non-superconducting metal 
YbCoIns have been grown, such that inversion symmetry is broken in the CeCoIns 
layers nzH]. This material is predicted to realize a nodal NCS with nondegenerate 
flat bands |44l [99]. Another intriguing compound is Cua;(PbSe) 5 (Bi 2 Se 3 ) 6 , where 
the doped topological insulator CuxBi 2 Se 3 is combined with the topologically triv¬ 
ial insulator PbSe m- Cu 2 ,Bi 2 Se 3 has attracted much attention as a candidate 
topological superconductor with a full (odd-parity) gap [1821 I183L I184j . Surpris¬ 
ingly, specific heat measurements appear to indicate that the superconducting gap 
in Cua;(PbSe) 5 (Bi 2 Se 3)6 has line nodes. Note that we have conservatively classified 
the pairing as spin singlet, although the upper critical field hints at a more exotic gap 
symmetry. 


7. Summary and outlook 


In this review we have presented an introduction to the topological properties of 
nodal superconductors. Although the topology is less robust than in fully-gapped 
superconductors, nodal superconductors are more common than unconventional fully- 
gapped superconductors, and so this theory will likely find much application. In 
addressing this subject, we have introduced two important examples: time-reversal- 
symmetry-breaking Weyl superconductors (Sec. 2.1), and time-reversal-invariant 
noncentrosymmetric superconductors (Sec. |2.2[), which can be regarded as prototypes 


of superconductors with topologically nontrivial point and line nodes, respectively. 

At the heart of the review is Sec. [^ and more specihcally the topological 
classification of nodal superconductors summarized in Table We have organized 
our review around this classification scheme, since: (i) It brings order to the growing 
zoo of unconventional superconductors; (ii) it gives guidance for the search and design 
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of new topological states; and (iii) it links the properties of the surface states to 
the bulk nodal structure and the symmetries of the order parameter. Point (iii) is 
particularly exciting, as it means that the pairing symmetry of an unconventional 
superconductor can be identified solely by means of surface measurements. Accurate 
modelling of the surface state spectrum is therefore an important aspect of the theory. 
This is discussed from the point of view of the popular quasiclassical approximation in 
Sec. 1^ While the theory agrees perfectly with the topological arguments, it also gives 
us access to physics beyond Table e.g., edge states without topological protection, 
and the physical attributes of all types of surface states. The latter is very useful 
input for designing experimental probes of the surface. 

In Sec. we briefly surveyed a number of materials for which there is strong 
evidence of a nodal superconducting state with nontrivial topology. Demonstrating 
the existence of the topologically protected edge states, however, will be a significant 
experimental challenge. As we have discussed in Sec. there have been many 
theoretical studies of the possible signatures of topological edge states, but only a few 
materials have been practically studied. While some systems are likely less suitable 
for such studies than others due to material-specific difficulties, the large number and 
diversity of candidate superconductors gives hope that direct detection of edge states 
is feasible. In this regards, proposals to realize a topological nodal superconducting 
state in heterostructures are very attractive, due to the fine experimental control 
over such systems. We expect this to be a major direction for both theoretical and 
experimental work. Another promising avenue for future research is to utilize the 
edge states of topological superconductors to design novel quantum devices, which 
might find technological applications in spin-eletronics and the development of new 
low-power electronics. A crucial first step towards the realization of such devices is 
the ability to control and manipulate the edge states using, for example, magnetic 
fields [8T1 II251 Il71j . electric fields |185[II86) . or mechanical strain. 

A significant omission in our discussion is the effect of interactions beyond their 
role in generating the superconducting gap. The BdG Hamiltonian, on which the 
topological classification is based, is a single-particle theory of the gapped state. 
Nevertheless, many candidate topological superconductors are strongly correlated 
materials, where interactions may exert a profound influence over the low temperature 
physics. The topological properties of interacting systems are of much interest, but 
remain poorly understood |187111881118911190j . A complete classification of interacting 
systems in one dimension has been obtained using matrix product states |191lll92lll93[ 
m, but such a classification is currently missing for higher dimensions. Although 
BdG theory generally gives a good description of unconventional superconductors 
(at least at a phenomenological level), interactions could still qualitatively alter the 
edge physics: For example, residual interactions may cause an instability of the 
surface states, in particular the flat bands with their divergent density of states. 
So far, attention has focused mainly on instabilities of the surface states of d-wave 
superconductors [1081II9511196L II97j , although the instability of flat bands of spinless 
superconductors has also been studied |I98) . 

To conclude, the discovery that electronic systems can have nontrivial topology 
has opened a new perspective on the fundamentals of condensed matter physics. Our 
review highlights the large changes wrought by this revolution in the subfield of nodal 
superconductivity, producing a new classification of this disparate class of materials 
by topological features of their gap structure. This does not supersede previous 
classifications, but rather complements them, and indeed it is already expected that 
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workers in this field pay attention to topological properties. We can only see this 
trend continuing, and that it will not be long before consideration of topological 
aspects will become ubiquitous and unremarkable in the study of unconventional 
superconductivity. 
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